arXiv:astro-ph/0211299vl 13 Nov 2002 


The dusty environment of quasars . Far-IR properties of Optical 

Quasars ^ 

Paola Andreani^ 

INAF - Osservatorio Astronomico Padova, Vicolo dell’Osservatorio 5, Padova, 1-35122, 

Italy, e-mail:andreani@pd. astro.it 

Stefano Cristiani^ 

INAF - Osservatorio Astronomico di Trieste, Via G.B. Tiepolo 11 - I 34131 Trieste, Italy, 

e-mail: cristiani@ts. astro.it 

Andrea Grazian^ 

Dipartimento di Astronomia, Universitd di Padova, Vicolo dell’Osservatorio 2, Padova, 

1-35122, Italy 

Fabio La Franca 

Dipartimento di Fisica, Universitd degli Studi ’’Roma Tre”, viale della Vasca Navale 84, 
I-OOI 46 Roma, Italy, e-mail:lafranea@fis.uniroma3.it 

Pippa Goldschmidt 

Astrophysics Group, Blackett Laboratory, Imperial College of Science Technology and 
Medicine, Prinee Consort Road, London SW7 2BZ, UK 

ABSTRACT 

We present the ISO far-IR photometry of a complete sub-sample of optically 
selected bright quasars belonging to two complete surveys selected through mul¬ 
ticolour (U,B,V,R)I) techniques. The ISOPHOT camera on board of the ISO 
Satellite was used to target these quasars at wavelengths of 7.3, 11.5, 60, 100 and 
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160 /xm. Almost two thirds of the objects were detected at least in one ISOPHOT 
band. The detection rate is independent of the source redshift, very likely due 
to the negative K-correction of the far-IR thermal emission. More than a half of 
the optically selected QSOs show signihcant emission between 4 and 100 fim in 
the quasar rest-frame. 

These fluxes have a very likely thermal origin, although in a few objects an 
additional contribution from a non-thermal component is plausible in the long 
wavelength bands. In a colour-colour diagram these objects span a wide range of 
properties from AGN-dominated to ULIRG-like. The far-IR composite spectrum 
of the quasar population presents a broad far-IR bump between 10 and 30 yum 
and a sharp drop at A >100 yUm in the quasar restframe. The amount of energy 
emitted in the far-IR, is on average a few times larger than that emitted in the 
blue and the ratio increases with the bolometric luminosity. Objects with 
fainter blue magnitudes have larger ratios between the far-IR (A > 60/um) fluxes 
and the blue band flux, which is attributed to extinction by dust around the 
central source. No relation between the blue absolute magnitude and the dust 
colour temperature is seen, suggesting that the dominant source of FIR energy 
could be linked to a concurrent starburst rather than to gravitational energy 
produced by the central engine. 

Subject headings: AGN - active galaxies: ISM, photometry - ISM: dust, extinc¬ 
tion, - Quasars: general - continuum - infrared: galaxies 


1. Introduction 

The properties of the circumnuclear environment are a key element to understand the 
formation and evolution of AGN. In particular, while these properties are well studied in 
the UV/optical spectral domain their mid-IR and far-IR characteristics are still a matter of 
debate. 

The sampling of the far-IR part of the Spectral Energy Distribution (SED) is important 
for determining not only the overall energy balance but also to discriminate the dominant 
emitting mechanism and to trace any evolutionary path from a dust-enshrouded phase to a 


^based on observations with ISO, an ESA project with instruments funded by ESA Member States (espe¬ 
cially the PI countries: France, Germany, the Netherlands and the United Kingdom) with the participation 
of ISAS and NASA 
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dust-free regime, corresponding to a ’’neat” optical phase (Sanders et ah, 1989). Possibly 
related is the investigation of the structure and anisotropy of the obscuring-reprocessing torus 
as predicted by unihcation models of AGN. In particular, knowledge of the spectral region 
where the dust emission peaks and starts to become optically thick is needed to constrain the 
nature of the torus and its fundamental geometrical and physical parameters (dimensions, 
dust temperature and optical depth) providing an orientation-independent parameter to test 
the unification hypothesis. 

IRAS observations have offered the hrst glimpse to the far-IR properties of optical 
quasars but did not have the required sensitivity to yield an unbiased view of the AGN class 
as a whole. In fact, only known optical quasars (PG quasars), with a very narrow range 
of fluxes and consequently a narrow range of luminosities at any given redshift, have been 
targetted and studied. The lack of any statistically well-dehned sample prevented from any 
inference on the general FIR properties of optically-seleeted quasars. Also snbsequent stndies, 
based on mm emission of IRAS QSOs and of high-z optical QSOs (Ghini et ah, 1989, Omont 
et al. 1996, Andreani et ah, 1999, Garilli et ah, 2000, Omont et ah, 2001, and references 
therein) and on ISO observations (Haas et ah, 1998, 2000; Polletta et ah, 2000; Oyabn et 
ah, 2001), could not address this issne and mainly focussed on the question of the emitting 
mechanism. 

The sampling of the shortest and longest parts of the IR spectrum via ISO observations is 
a key test to disentangle the emitting mechanism. However a typical conclusion of ISO studies 
is that it is not easy to disentangle pnre AGN emission from a likely common FIR-emitting 
star-formation component. There is not a general agreement among different observations 
(Haas et ah, 1998; 2000; Polletta et ah, 2000) and interpreting models (e.g. Granato & 
Danese, 1994; Rowan-Robinson, 1995). 

In this paper we report ISO observations of a complete optically selected sample of bright 
(15< mB <17) qnasars. The original aim of our programme was to use these data to infer 
any evolutionary link from ’’pure” AGN objects to eventual transition objects linking AGN 
to ULIRGs. The investigation of the complete sample would have allowed us to perform a 
number of statistical tests (for instance the evaluation of the bivariate Inminosity fnnction 
at optical-IR wavelengths) to stndy how the optical-infrared relation for qnasars changes as 
a fnnction of both luminosity and redshift. Because of the reduced sensitivity of ISOPHOT 
the number of observed objects is greatly reduced from the original one and the number of 
detection is also consequently too small to allow reliable statistics. Neverthless, a number of 
different properties of the far-IR emission of optical quasars can be inferred from these data 
and help nnderstanding the expected behaviour of the type-1 AGN population in surveys 
from future space mission - such as SIRT-F and HSO. 
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2. The Sample 

A statistically well defined sample of bright quasars was selected from two complete 
optical multicolour (U,B,V,R,I) quasar surveys: the Edinburgh Bright Quasar Survey (Gold¬ 
schmidt et ah, 1992) and the ESO Key-Programme ”A Homogeneous Bright Quasar Sur¬ 
vey” (Cristiani et ah 1995) in a total sky area of 888 sq. deg. The optical fluxes are in the 
range 15 < B < 17, partially overlapping those of the PG survey, but with a much better 
photometric accuracy {as ~0.1mag). The objects finally observed with ISO are listed in 
Table 1, where names (column 1), J2000 coordinates (2-3), redshifts (4), Johnson Consins 
U,B,V,R,I magnitudes (columns 5-9), and the absolute B magnitudes are reported. The 
apparent B-magnitude, me, were taken from the Homogeneous Bright QSO survey (HBQS) 
(Gristiani et al. 1995) whereas Bj apparent magnitudes were derived from the Digitized Sky 
Survey (DSS) ^ calibrated plates (Grazian et ah, 2001). K-corrections were inferred from a 
QSO composite spectrum (Gristiani and Vio 1990). The apparent magnitudes are corrected 
for Galactic absorption using the reddening maps derived from HI and Galaxy Gounts by 
Burstein & Heiles (1982) and the luminosity distances are derived according to the adopted 
Gosmology (Ho=50 km/s/Mpc, A = 0, D = 1). 


3. The Observations 
3.1. The ISOPHOT Observations 

The ISOPHOT data presented here result from the merging of two distinct Open Time 
programmes PANDREAN and PGOLDSGH. 

The observed sample consists of 34 quasars of which 21 were observed at 11.5, 60 and 160 
/rm and 13 at 7.5, 25, 60 and 100 /im. The log of the ISOPHOT Observations is in Table 2, 

coordinates nsed for pointing are those listed in Table 1. 

The data were reduced with the ISOPHOT Interactive Analysis Tool (PIA Version V8.2, 
Gabriel et ah, 1998). All data were corrected for electronics nonlinearity and cosmic particle 
glitches. Ghopped data were then averaged: all the points corresponding to one chopper 
position were averaged and fitted with a polynomial curve to estimate their time behaviour. 
If the detector shows stability, i.e. the fitted curve converges to a constant value, this latter is 
used as the signal related to that chopper position. In some cases, when the integration time 


^The online Digitized Sky Surveys (http://archive.eso.org/dss/dss) server at the ESO/ST-ECF Archive 
provides access to the CD-ROM set produced by the Space Telescope Science Institute through its Guide 
Star Survey group. 
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on source was too short, for the ISOPHOT detectors to reach the stability, the signal value 
was inferred by extrapolating the time-dependent signal curve to an asimptotic constant 
value. Data are further corrected for glitches, for orbital dependent dark current and for the 
signal dependence on the ramp integration time. 

The same reduction procedure was applied to the associated measurements of the thermal 
Fine Calibration Source (FCSl), which were used to estimate fluxes. 

For raster maps additional softwares were used to extract the signal, kindly provided by 
Use van Bemmel and Martin Haas and are described in van Bemmel et ah (2000) and Haas 
et ah (2000). 

Table 3 lists the ISOPHOT fluxes at 7, 12, 25, 60, 100 and 160 fim (columns 2-7). 

30 % of the objects have fluxes larger than 3 cr at 60 ^um, 5 out of 7 were detected with raster 
maps while only 5 out of 27 with chopper mode. No objects were detected at 7 and 12 jim 
(only in chopping mode) and only 4 at 25 and 100 /xm. At 160 /xm 6 objects were detected 
(3 detections obtained from the 7 raster scannings). The raster scanning observing mode 
shows the highest sensitivity. However at 160 /xm the sensitivity could be reduced in regions 
where the Galactic background due to cirrus emission is high. We checked detections against 
cirrus emission level and the results are shown in the Appendix. 

Figure 1 provides a general view of all the ISOPHOT fluxes collected in the present 
work. The left panels show them as a function of redshift, the right panels as a function of 
the B-band flux. MIR is the observed wavelength range 7-25 /xm and FIR 100- 160 /xm. 
The lack of detections, shown in Figure 1, can be ascribed not only to the reduced sensitivity 
in the chopper mode but also to the nature of the FIR emission in optical quasars, showing a 
dehcit in the short (7 and 12 ^m) and long (160 /xm) observed wavelengths emission, because 
of the characteristic high temperature of the dust grains heated by the AGN. 

Although many of the values reported in Figure 1 are upper limits, it emerges an absence 
of dependence of the detection rate with redshift: the detections are randomly spread all 
over the redshift range. 10 detected objects have redshifts larger than the mean redshift 
{{z) = 0.72, median redshift is 0.7) of the sample: 2 : > (z), while the other 10 detections 
have < (z). Similarly there is little correlation between the B apparent luminosity and 
the FIR fluxes, with only a slight trend of larger FIR fluxes for smaller B flux. The most 
B-luminous objects are not correspondingly FIR-luminous and there is no scaling of the 
FIR fluxes according to the B-band magnitude as it was assumed computing the integration 
time (the smaller the magnitude the shorter the integration time). The latter was computed 
aiming at a uniform sensitivity in terms of aoiR = logFbiue — log Fir ~ —2. 

Only four objects are detected in more than one IR band: 0144-39, 1404-1-09, 1415-1-00 
and 1415-00. Individual notes are given in Section 4. 
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redshift Fg (mJy) 


Fig. 1.— Logarithm of the ISOPHOT fluxes in mJy as a function of redshift (left panels) 
and of the B fluxes in mJy (right panels). The lower panels report the fluxes at 7 and 12 
/im (circles) and 25 /rm (diamonds), the upper panels the fluxes at 100 /rm (triangles) and 
160 /im (pentagons). Filled symbols correspond to detections, open to upper limits. No 
clear dependence is seen: upper limits and detections are equally distributed (see text). A 
slight trend - lower IR fluxes at large B fluxes - exists for those objects detected at 60 and 
160 /im. 
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3.2. Additional data 

3.2.1. IRAS SCANPI data 

The IRAS data listed in table 3 are the result of co-added fluxes provided by IPAC 
SCANPI (Scan Processing and Integration Tool) program. This procedure performs a one¬ 
dimensional coaddition of all the IRAS survey data on the source. The sensitivity is com¬ 
parable to that achieved by the FSC (Faint Source Catalog) for point sources (see the IPAC 
manual for details, http://ipac.caltech.edu/). 

Tentative (3cr) detections among IRAS data are found for 0120-28, 0144-39, 1252-1-02, 
1321+28, 1326-05, 1351+01, 1355+02, 1415+00, 1415-00, 2313-30, 2357-35 and are shown 
in boldface in table 3. 

Some objects in table 3 show discrepant 60+m fluxes between the ISO and IRAS 
measurements. Some authors (e.g. Alton et ah, 1998; Spinoglio et ah, 2002) have 
argued that ISO fluxes are affected by a calibration error as large as 30% and 
found ISO fluxes in excess with respect the corresponding IRAS ones. However 
in some cases (1415-00 and 1415+00), even taking into account an additional 
30 % error ISO values lie well above the IRAS ones. As shown in Appendix A, 
it is also not straightforward to attribute this discrepancy to additional cirrus 
emission. This contribution is expected to be larger in IRAS measurements be¬ 
cause of the wider IRAS beams. This discrepancy may be due to variability but, 
at present, this issue cannot be settled and could only be addressed by regular 
radio observations. 


3.2.2. near-IR Photometry 

Counterparts within 4 " of the QSO listed in Table 1 were detected using the 
2MASS Second Incremental Data Release (http://www.ipac.caltech.edu/2mass/) 
and the corresponding J (1.25+m), H (1.65+m) and (2.17/xm) photometry of 
the 19 (out of 34) objects found is listed in table 1 (see also Barkhouse &; Hall 
2001). All objects with near-IR counterparts show rising near-IR spectra with 
magnitudes decreasing from the J to the K-band. 



3.2.3. radio data 


Radio data are used in this work to check whether strong non-thermal emission showing 
up in the radio domain would possibly affect the far-IR ffuxes. The data were taken from the 
NVSS (NRAO VIA Sky Survey) survey (Condon et al. 1998) and from the FIRST survey 
(Becker et al.l994) and are listed in table 1 for the 21 objects with sky positions covered 
by these surveys (at 6 > —40°). 8 objects out of 21 observed in these surveys have 
a logarithmic ratio between the 1.4 GHz and the B-band luminosities, 
larger than 1, i.e. are radio-loud objects. 


4. Results 

4.1. Notes on individual objects 

• 0051-39: has only one detection at 160 /rm. The low value of the 60 iim upper limit 
suggests that the 160 /rm flux could be eventually affected by either a galactic cirrus 
or radio non-thermal emission. The lack of radio data on this object does not allow us 
to check the latter hypothesis. 

• 0059-30: has only one detection at 160 /rm. This object is radio-loud and therefore 
this measurement could be affected by a non-thermal component but its flux is still 
compatible with the emission of a putative host galaxy at this wavelength. 

• 0129-40: has a very clean 60 ^m detection, obtained with observations in raster mode, 
confirming the higher sensitivity of this option. 

• 0144-39: detected both at 60 and 160 ^m through raster pointings. The 60 ^m IRAS 
flux agrees with the ISO detection, while that at 100 ^m is also in agreement with that 
at 160 nm for a thermal spectrum varying with wavelength as 

• 1252-1-02: detected at 160 /im and marginally at 60 iim. IRAS found ffuxes at 12 and 
25 fim. It is very difficult to combine these four measurements in a unique thermal 
component. They likely belong to different thermal components. 

• 1355-1-02: detected at 60 iim and possibily at 100 frm by IRAS with an upper limit at 
160 /xm. Its spectrum can be reproduced by the emission of the circumnuclear torus. 

• 1404-1-09: has detections at 25, 60 and 100 frm. Since it is a radio-loud these measure¬ 
ments could be due to the combination of two different emission mechanisms, thermal 
and non-thermal (see also (Polletta et al. 2000)) 
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A(/xm) 


Fig. 2.— Spectral Energy Distributions (SEDs) as vL{v) normalized at 60 /im of the objects 
reported in tables l-i-3 having at least one ISOPHOT detection. Open squares refer to optical 
and near-IR photometry, hlled circles to ISOPHOT data, asterisks to IRAS data. The solid 
line corresponds to one of the models by Granato & Danese (1994) predicting the emission 
of dust distributed in the surrounding medium of the central AGN (see text). The dotted 
line is the spectrum of M82 shown here as a prototype starburst. 
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Fig. 4.— Logarithm of the ISOPHOT flux ratios versus redshift (left panels) and B apparent 
magnitude (right panels). The lower panels report the 7,12 and 25^m (MIR) fluxes, the 
upper panels the 100 and 160 fim (FIR) fluxes. Symbols are as in Figure 1. A weak redshift- 
dependence is seen in the left panels for long wavelength fluxes, while the ratio does indeed 
show a dependence on me. In the right panels upper limits and detections are equally 
distributed. The solid line in the upper left panel shows the expected behaviour of the ratios 
in presence of extinction. 
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• 1406-03, 1415-1-00 and 1415-00 have all large fluxes at 60 and 100 /rm. The latter two 
were detected also at 25 jim. We cannot exclude that the large 60 /xm fluxes are affected 
by instrumental artifact. Galactic cirrus can contribute to the 100 /im emission. The 
100 IRAS maps in these regions show a quite high background of around 80 MJy/sr 
with a variation from pixel to pixel around 10 %. This means that 3±0.3 Jy fall on 
the ClOO pixel. An imperfect beam switching could then mimic a source emission. 
However, we are not aware at present of such an effect in the ISOPHOT data. 

• 2313-30: the only object observed both in raster and chopping modes but detected 
only with the former instrumental set-up at 60 /im. The chopper data are noiser but 
consistent with the raster detection. 

• 2321-29 and 2357-35 have a detection only at 160 /xm. Without any other FIR mea¬ 
surements it is tough to infer the nature of this flux. It could be ascribed to a cool 
component of the host galaxy as, for instance, seen in radio galaxies (van Bemmel 
et ah, 2000). 


4.2. Spectral Energy Distribution 

Figures 2 and 3 show the spectral energy distribution in the object’s rest-frame of the 
18 sources in the sample with FIR detections. Spectra are normalized to the 60 /xm flux. 
Open squares correspond to the UBVRI photometry and to JHK data, filled circles to ISO 
measurements (down-arrows are upper limits), stars to IRAS points, which are in general 
upper limits but for 0105-26, 0120-28, 0144-39, 1252+02, 1351+01, 1355+02. 

Although the spread in the SED appears large, mainly due to the presence of censored 
data, some common features can be identified: (a) all quasars with good near-IR data show 
a dip in the SED around A = 1 /xm (b) all objects (but 0051-39) with 160 /xm data show a 
drop of the far-IR SED beyond 100 /xm. 

The comparison between the optical and far-IR spectrum must take into account possible 
variations due to the intrinsic variability of these sources. The effect of the variability 
is difficult to address in detail, as it significantly depends on the time lag between the 
various flux measurements, the quasar absolute luminosity and increases with decreasing 
wavelength. From the analysis of the structure function (Cristiani et al. 1996) we should 
expect an average uncertainty due to variability of 0.2 magnitudes for quasars having a 
typical absolute magnitude ~-26 (0.1 for ~-30 and 0.3 for ~-23). Such variations 
of the flux values are still within the size of the photometric points shown in Figures 2 and 
3. 
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The thick solid curves in Figures 2 and 3 represent the predictions of the Granato and 
Danese’s model (1994), which computes the expected emission of a circumnuclear torus-like 
dust distribution around the central nuclear source. The solution of the radiative transfer 
equation is required since in such tori the dust emission is self-absorbed even in the near- 
and mid-IR. 

Dust is illuminated by the central source and reaches an equilibrium temperature which 
is a function of the intensity of the radiation held and of the dimensions and chemical 
composition of the assumed grain components. In the case shown here the chosen geometry 
for the dust distribution has an axial symmetry and a scaleheight along the vertical axis 
increasing linearly with distance (‘hared disc’). The central source is either unobscured or 
only partially obscured from view. This model naturally explains the SED dip around 1 ^m: 
dust too close to the central engine attains too high temperature and sublimates. 

All the curves shown in Figures 2 and 3 correspond to the same ’face-on’ model, being 
all these UV-excess objects with the central source un-extincted by the surrounding torus 
material. The diherence in the amount of IR huxes can be accounted for by varying the dust 
distribution size, i.e. the dimension of the torus. 

We have neither attempted to hnd the best-ht curves, because of the lack of sufficient con¬ 
straints from the data, nor we claim that the model shown is the ’best’ to explain these 
observations. Other models would be equally good (see e.g. Rowan-Robinson, 1995). 
However, the comparison between the predicted curves and the photometric points provides 
useful insights: a simple torus-like dust distribution around the central source is a too naive 
picture. In some cases the required dimension of the dust distribution is wide and dust could 
be broadly spread in a putative host galaxy disc. 

To show how an additional galaxy disc emission modihes the SED, the optical-mm spectrum 
of M82, a typical local starburst, is also plotted in Figures 2 and 3, as a dotted line. 


4.3. The colour plots 

Figure 4 reports the colours, as a function of the B-band apparent 

magnitude (right panels) and of the source redshift (left panels). 

The increase of the ratios with the B-magnitude can be due to several effects: (a) an obser¬ 
vational bias, (b) cirrus emission or (c) extinction. 

(a) The upper left corner of every hgure panel would contain objects with large far-IR 
emission and bright apparent blue magnitude and therefore easily observable. Objects in the 
lower right corner could be missed because of a lack of sensitivity of far-IR observations. The 



14 



15.5 16 16.5 17 17.5 

m 

I I ig 


Fig. 5.— Flux ratios between the K-band and the B-band, expressed as hib — hik, against 
the B apparent magnitude for those objects with near-IR photometry reported in Table 1. 

The solid line corresponds to the estimated linear regression line, the standard 
deviation, 0.575, is large and the correlation coefficient is 0.33. The hypothesis 
that does not vary with me can only be rejected at a 13% level. 





F(60 /xm)/F(B) 


Fig. 6.— The ratio between the 100 and 160 /im fluxes and the B-band flux is shown as a 
function of the ratio between the 60 iim and the B flux. Open squares correspond to the 
100 /xm data while filled pentagons to the 160 fim ones. The straight line corresponds to 
a one-to-one relationship between the two colours. The outlier at ~ 600 and ~ 50 
corresponds to 0051-39, which, as shown in Figure 2, has an extraordinarly large 160/xm 
flux. 
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upper limits distribution is quite uniform with respect to m^, i.e. there is no accumulation of 
upper limits at the faint end. Sensitivity limits therefore do not seem to affect dramatically 
the results. 

If the variation of the ratios in Figure 4 were only linked to the QSO spectrum, we would 
expect to see this ratio decreasing with redshift: the B-flux increases with redshift because 
of the negative k-correction and for a constant FIR flux would be smaller at higher 
redshifts ®. The observations show the opposite trend. 

(b) The presence of Galactic cirrus emission falling in the ISOPHOT beam and affecting 
the far-IR fluxes has been checked using independent measnrements at these wavelengths. 
We used the IRAS 100 pm maps in regions surronnding the observed objects (see Appendix 
for detail). Although in three cases the far-IR fluxes can be really enhanced because of a 
cirrus contribution (see Appendix) there is not a one-to-one correspondence between the 
ratios and the 100 pm rms values. 

(c) The observed behaviour could be then ascribed to an extinction along the line 
of sight. The left npper panel shows a small increase of with redshift. The solid 
line corresponds to the expected change of this ratio when the emitted UV-optical flux 
(F(A) oc A“° '^) varies, because of the extinction A(A) and k-correction, k{z), as 

Fb = F(A,z)A;(z) ■ (1) 

We compute the ratio for blue fluxes, Fb, affected by an extinction of A(A) = Es{B — 
V) * Qexi(A). We adopted as extinction law that provided by Calzetti (2001) for starburst 
galaxies. 

Neglecting k-correction of the far-IR fluxes and correcting B-fluxes as in eq.l, more than two 
third of the ratios, in the upper right panel of Fig.4 would be shifted to an average 

valne of 100. This conhrms that the gross behaviour of is dictated by extinction. 

Onr interpretation is therefore that fainter objects in the optical are also redder. The 
redshift-dependence of IR to blue flux ratios were also investigated by Alexander et al. 
(2001) for a sample of ELAIS AGN with 15 pm fluxes. These authors find a ratio 
ranging between 2.5 and 100 independent of the source redshift. Haas et al. (2000) as well 
do not detect any variation of both and with redshift. 

To further check the presence of dust extinction along the line of sight throngh the sources, we 


®k-corrections of FIR fluxes are neglected here since, within the redshift range of our objects, the observed 
FIR fluxes sample the Wien region of the thermal spectrum and the resulting variation of the FIR fluxes in 
the restframe and observed frame is small. 
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plot in Figure 5 the ratio between the K-band and the B-band fluxes, expressed as uib — niK 
against the B magnitude. The straight line in Fignre 5 represents to the estimated 
linear regression line. In spite of the large spread (standard deviation is 0.575), 
there is a marginal correlation (the correlation coefficient is 0.33) and a slight 
probability that at larger values of me — mx correspond larger m^ in accordance 
with the ratios plotted in Figure 4. However, the hypothesis that does not 
vary with me can only be rejected at a 13% level. 


^.3.1. Colour-colour plot 

In Figure 6, 100 and 160 pm fluxes, normalized to the B-band flux, are shown 
against the 60 pm flux, The two colours are tightly related: an increase of the long 
wavelength flux with respect to the B-band flux corresponds to an equivalent enhancement 
of the 60 pm flux with respect to the optical one. The two far-IR fluxes are therefore related 
and very likely belong to the same thermal component. As a comparison a straight line 
showing a one-to-one relationship is drawn in Figure 6. The ratio can be then used 

as an estimation of the dust colour temperature for that thermal component. 


4.4. Dust Temperature 

Customarily the spectral indices, a, for a oc z/" SED are defined as 

tt(z/i,z/ 2 ) = log^^^ - log— (2) 

and simply related to colours. Figure 7 reports Q;(60^m, MIR) versus q:(FIR, 60/rm). Ac¬ 
cording to our previous finding (§4.3) this latter is inversely proportional to the warm dust 
temperature, T^d- Although many points in the diagram are upper/lower limits there seems 
to be a trend: the relative emission at shorter wavelengths (7, 12 and 25 pm) with respect 
to the 60^ emission decreases as T^d increases. This colour-colour diagram was widely used 
in the past as a tool to detect and discriminate between different types of activity in the 
nuclear and circumnuclear regions of galaxies (see e.g., Canalizo & Stockton (2001) and 
references therein). Different kinds of objects (QSOs/Seyferts, starbursts and ULIRGS) oc¬ 
cupy distinct loci in such a diagram: objects for which FIR fluxes are dominated by dust 
reradiation are located in the lower right corner (in our convention at lower T^d and lower 
Q;(60^m, MIR)), while objects with strong non-thermal emission, such as optically selected 
QSOs, are in the upper left corner. Indeed PG quasars observed by IRAS are preferentially 
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found at —2 < Q;(60^m, MIR) < 1 and —2 < q:(FIR, 60/im) < 0 partially overlapping the 
region of the diagram where transition objects are found (Canalizo & Stockton, 2001). To 
show that in detail we plot in Figure 7 the location of the QSOs observed by IRAS at redshift 
lower than 2 and studied by Andreani, Franceschini & Granato et ah (1999). 

In Figure 8 the ratio is plotted against the absolute B magnitude. With all the 

limits of this analysis due to the poor information on the FIR fluxes, it appears that no 
statistically signihcant correlation is present between the B luminosity and the warm dust 
colour temperature, T^d- 

If Mb is linked to the nuclear emission, while is only a dust property, it seems that no 

strong relation exists between the energy emitted by the nuclear source and that emitted in 
the FIR. Similar hndings were already reported by Andreani, Franceschini & Granato (1999) 
for an inhomogeneous sample of optical quasars with IRAS and mm-wavelengths fluxes. For 
their sample no relation was found between the dust mass and the absolute B-magnitude 
Mb. McMahon et ah (1999) and Omont et al. (2001) do not detect any dependence of the 
sub-mm fluxes on the optical luminosity in high redshift quasars. All these results could 
be htted into a picture in which the emission from the nuclear source is not tightly related 
to the physics (the dust properties: temperature, mass and luminosity) of the surrounding 
medium and argues in favour of the hypothesis that most of the FIR luminosity is linked to 
a concurrent starburst and trace activity in the host galaxy (see also Glements, 2000). This 
result is also related to dynamical studies aiming at the search for BHs in galaxy centres. 
These investigations show that although there is a strong one-to-one correspondence between 
the central BH mass and bulge luminosity no detectable correlation is found between the 
BH mass and the luminosity of galaxy discs (see e.g. Kormendy & Gebhardt,2001). 


4.5. Luminosity plots 

Figure 9 shows the ratio between the far-IR, LpiR, and the B-band, Lb, luminosities ver¬ 
sus the bolometric luminosity, Lboi- On the y-axis we plot i-6- fhe ratio between the 

far-IR luminosity LpiR, as the sum of the different far-IR data, and the blue-band luminosity 
computed from the absolute B magnitude: 10“0 '^'^b+28.523 Qj^ x-axis the bolometric 
luminosity is shown: 


Lboi = 4:7rdl 


du /(z/) 


( 3 ) 


where di is the luminosity distance. Lboi is computed from UV to FIR. Down-arrows corre¬ 
spond to objects undetected in all bands, triangles to those with one detection, hlled circles 



19 



-2 0 2 


a(FIR,60) 

Fig. 7.— The ratio between 60^m flux and MIR fluxes is shown as a function of the ratio 
between FIR and 60 /rm fluxes. Ratios are defined by equation 4.4. Asterisks correspond to 
the 25 /iui data, filled circles to 12 /rni, open circles to 7 ones. Open squares correspond 
to QSOs at z < 2 reported in Andreani, Franceschini & Granato (1999). Distinct areas 
can be identified from AGN dominated (in the upper left corner ) to Starburst dominated 
objects (in the lower right corner), according to the expected dominant emitting mechanism 
in the object SED from non-thermal to thermal. 





-24 -26 -28 -30 


Mb 

Fig. 8.— The ratio , which is is proportional to the dust temperature, T„d, as a 

function of the absolute B magnitude, Mb- Asterisks refer to the 100 pm data, hlled circles 
to the 160 pm data, © symbols to ratios with upper limits on both fluxes. 
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0 0.5 1 1.5 2 


redshift 

Fig. 9.— The ratio between the far-IR, Lfir, and the B-band, Lb, luminosities against 
the bolometric luminosity, Lboi (upper panel) and against the source redshift (lower panel) 
is shown. Triangles correspond to those objects with only one ISOPHOT detection, hlled 
circles to those with two or three detections, downarrows to undetected objects. 
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to those detected at least in two bands. Lfir could well be underestimated for those objects 
with only one or two detected fluxes. Sampling the Wien region of a thermal spectrum with 
a characteristic temperature of 40-100 K at only one or two frequencies means to underesti¬ 
mate the total luminosity of a factor between 2 and 5 with respect to the integral over the 
whole frequency range. 

Only hve objects have clearly a ratio lower than 1, 14 detected ones have > 1 
of which 7 have FIR luminosity more than one order of magnitude higher than the blue 
one. There is a clear trend of increasing with increasing bolometric luminosity. This 
means that the FIR luminosity becomes increasingly important as the bolometric luminos¬ 
ity increases. A similar behaviour is seen in nearby ultraluminous IR Galaxies (ULIRGs) 
(Sanders and Mirabel, 1996) and has been associated to the increased central concentration 
of molecular gas. Because of the high extinction the relative role of nuclear starburst and 
AGN activity is hard to disentangle in ULIRGs but there seems to be an enhancement of 
AGN contribution to the overall energy budget as the bolometric luminosity increases. Fur¬ 
thermore, in the Unihed Scheme of AGN an enhancement of the ratio is expected, since 
the reprocessed emission by the dusty torus emerges isotropically in the FIR. The colour axis 

should be related to the viewing angle to the torus with more inclined objects having 
larger therefore lying on the right-hand side of the diagram. 

The lower panel in Figure 9 shows the behaviour of the ratio against the 
source redshift: with the present data it is not possible to detect any evolution 
of the ratio with redshift. 


5. The quasar composite spectrum 

The objects presented in this work are randomly drawn from a complete sample of 
optically selected quasars (see §2) and can be considered representative of the entire quasar 
population. If we assume that the radio-quiet quasars can be modeled as a homogeneous 
population, a spectrum can be built via observations at hxed frequencies of targets at different 
redshifts. This analysis exploits the whole spectral information from the optical to the FIR 
and helps outlining possible general features of the emitting mechanisms and the physics of 
the QSO environment. The spectrum is built by dividing the available wavelength range in 
bins. In each bin an average of the fluxes at that wavelength is computed. The errorbars 
shown are those related to the averages. Uncertainties on the flux scale lie well within the 
large errorbars. 

Figure 10 shows the average spectrum in the QSO restframe. The spectrum is nor- 
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malized at Ab = 0.44 yum. Averages are computed in each wavelength bin with the survival 
statistics which take into account the censored data (Feigelson 1990). Detections with a 
suspected contribution from IR-cirrus were included as upper limits. A well dehned IR com¬ 
ponent peaking at 10-30 //m and dropping steeply above 100 ^m is evident. The spectral 
shape of the far-IR component is characteristic of a thermal origin. It can be speculated 
whether it is due to starburst emission in the host galaxy and/or to a dusty torus around 
the central source. In order to derive constraints to the available parameters on the dust 
distribution, we can compare the data with available models. Figure 10 reports three differ¬ 
ent spectral energy distributions predicted by the Granato & Danese ’s model (1994) (see 
§4.2). We restrict ourselves to models which allow UV-optical photons to escape from the 
nuclear source since the quasars were selected on the basis of their UV-excess. The dust dis¬ 
tribution around the central source should be compact because of the sharp cut-off at long 
wavelengths of the FIR bump. A larger bump would imply larger temperature distribution 
and very likely larger spatial dimension of the dust distribution. 

The solid line corresponds to face-on (with a ’naked’ nuclear source) compact (i.e. the 
ratio between the inner and outer radius of the torus ^^=800) conhguration of the torus 
and an optical depth r at 0.3 jim of 30; the dotted line to a 45° inclined torus with larger 
dimension (/^=1000) and same r, the dashed line to a 45° inclined torus (/^=1500) with 
higher absorption, r=60. The curve in better agreement with the available photometric 
points is that corresponding to the larger optical depth (dashed line) and therefore the 
strongest dust self-absorption around 10-30 /im. 


6. Conclusions 

A sample of optically selected quasars was observed with ISOPHOT between 7 and 160 
yum, with the following results: 

• The sample shows homogeneity in mid and far-IR properties: the mid and far-IR 
colours depend neither on the redshift nor on the absolute B-magnitude and there is 
no correlation between the blue magnitude and the ISOPHOT fluxes. 

• The detection rate is not a function of redshift ( 10 detected objects have < 0.7, the 
other 10 have 2 : > 0.7, where 0.7 is the median redshift of the redshift distribution) 
consistently with the negative K-correction implied by a thermal spectrum. 

• The long wavelengths (A=60,100 and 160 yum) fluxes are tightly related and likely arise 
from the same thermal component. 




og A(/[xm) restframe 


Fig. 10.— The average spectrum in the QSO restframe. The spectrum is normalized at 
Ab = 0.44 //m. Filled circles correspond to the optical and near-IR data, open triangles to 
the ISOPHOT measurements discussed in this paper. The lines going through the data show 
the predictions of the models by Granato and Danese (1994) which are able to £t the QSO 
SED (see §4.2). The solid line corresponds to a face-on compact conhguration of the torus 
and absorption at 30/xm of r=30, dotted line to a 45° inclined torus with larger dimension 
and same r, the dashed line to a 45° inclined torus with higher absorption, r=60. 
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• Objects with fainter blue magnitudes, me, have larger ratios between the FIR (A > 

60/im) fluxes and the blue band flux, We ascribe this behaviour to a larger 

extinction along the path to the central source. 

• The colour-colour diagram a(60/im, MIR) versus a(FIR, 60/im) shows that the objects 
of the present sample span a wide range of properties from AGN-dominated ones (low 
q;(FIR, 60/im), i.e. large dust temperature) to starburst-dominated ones (low and 
large ^). 

• No statistically signihcant correlation is detected between the absolute B-magnitude, 
Mb and the warm dust temperature, T^d- If Mb mainly reflects the nuclear emission 
and the FIR fluxes the emission of the surrounding dust, it seems that no relation 
exists between the energy emitted by the nuclear source and that by dust. The latter 
could rather be a tracer of the starburst activity in the host galaxy. 

• Even in an optically selected sample the amount of energy emitted in the far-IR, Lfir, 
is on average a few times larger than that emitted in the blue. Lb. The ratio does 
not show any clear variation with redshift but increases with the bolometric luminosity, 
Lboio- This behaviour is similar to that seen in local ULIRGs, possibly because of an 
increased central concentration of molecular gas. In the Unihed Scheme of AGN an 
enhancement of the ratio is expected in the more inclined objects for which the 
viewing angle goes through the torus, i.e. through a more extincted line of sight. 

• The QSO Gomposite Spectrum, i.e. the average spectrum built in the QSO restframe, 
shows a broad far-IR bump at 10-30 /xm which could be due either to a starburst 
emission in the host galaxy or to a dusty torus around the central source as predicted 
in the unihed scheme of AGN. 

Restricting our analysis to only the QSO SED one could argue that the best picture 
htting the data is that of a circumnuclear dusty torus. However other arguments 
presented also in this paper (the colour-colour diagram, the lack of relation between 
the power of the nuclear source and the dust temperature and the increase in the ratio 
see Figures 7, 8, 9) suggest that a contribution from a starburst from the host 
galaxy is quite unavoidable. 


This publication makes use of data products from the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Genter/Galifornia Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. 
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Telescope was operated by the Royal Observatory Edinburgh, with funding from the UK 
Science and Engineering Research Council, until 1988 June, and thereafter by the Anglo- 
Australian Observatory. Original plate material is copyright the Royal Observatory Edin¬ 
burgh and the Anglo-Australian Observatory. The plates were processed into the present 
compressed digital form with their permission. The Digitized Sky Survey was produced at 
the Space Telescope Science Institute under US Government grant NAG W-2166. 

One of the authors (PA) warmly thanks Use van Bemmel for her help during the data anal¬ 
ysis and for her permission of using her software. PA acknowledges the fundamental help 
of the ISOPHOT Data Center in Heidelberg in the data reduction and MPE for hospitality. 
Alexander von Humboldt Foundation for support is acknowledged. Part of this work was 
funded by the Italian Space Agency (ASI) under contract ARS-98-226, ASI-I-R-105-2000. 


A. Appendix: Detections against cirrus level 

Figure 11 reports the detected fluxes and upper limits at 100 and 160 pm against the 
r.m.s. values of the cirrus emission as estimated from the 100 pm IRAS maps in the sky 
positions corresponding to the source location. Upper limits are shown by downarrows. 
No correlation between the two quantities is apparent, with the possible exception of the 
three sources 0051-39, 1404-1-09, 1406-03. For these three cases the large observed flux could 
be affected by cirrus emission as discussed in section 4.1. However, 0051-39, 1404-1-09 are 
radio-loud objects and the tail of the radio emission could also contribute to the FIR fluxes. 
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at 100 ^m as estimated from the IRAS maps towards the sources observed by ISO. Upper 
limits are indicated as downarrows. The three higher values of the ISO fluxes correspond to 
0051-30, 1404-1-09 and 1406-03, the first two are also radio-loud objects. 
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Table 1. Optical and near-IR Photometry of Quasars 


Name 

R.A.(J2000) 

dec(J2000) 

z 

me 

mu 

mv 

mR 

mi 

mj 

mn 

niK 

Mb 

0034-70 

0 34 05.0 

-70 25 54 

0.363 

15.30 


15.50 

17.70 


16.10 

15.61 

14.59 

-26.27 

0049-29 

0 51 38.1 

-29 23 13 

0.308 

15.84 

16.66 

14.81 

13.96 

13.59 

14.99 

14.37 

14.89 

-25.40 

0051-39 

0 54 09.5 

-39 16 52 

0.224 

16.93 

16.92 

17.30 

16.62 

16.51 




-23.63 

0059-30 

1 02 00.7 

-30 18 28 

1.033 

16.92 

16.41 

17.01 

16.64 


16.11 

15.77 

15.26 

-26.98 

0105-26 

1 05 58.9 

-26 06 22 

0.776 

16.60 


16.90 

15.60 





-26.65 

0120-28 

1 22 36.6 

-28 43 20 

0.434 

16.05 

15.48 

15.71 

15.84 

15.91 

15.48 

14.68 

13.84 

-25.91 

0129-40 

1 31 43.9 

-40 36 54 

1.371 

15.75 

15.06 

15.72 

15.13 

15.16 

14.68 

13.89 

13.71 

-28.81 

0144-39 

1 46 12.4 

-39 23 06 

0.244 

15.88 

15.57 

15.87 

15.56 

15.07 




-24.88 

0203-08 

2 03 22.5 

-8 43 37 

0.770 

16.20 

16.50 

16.90 


15.84 

15.48 

15.12 

-27.03 

<24.5 

1245-03 

12 47 35.0 

-3 50 09 

0.379 

15.93 

15.04 

16.07 



15.41 

14.80 

13.77 

-25.73 

1249-02 

12 51 51.4 

-2 23 36 

1.184 

17.00 

16.45 

17.10 



15.80 

15.30 

14.84 

-27.19 

1249-06 

12 51 56.4 

-7 05 02 

1.187 

16.58 

15.22 




15.52 

14.99 

14.55 

-27.61 

1252+02 

12 55 19.6 

1 44 13 

0.345 

15.52 

14.13 

15.48 






-25.95 

1321-05 

13 24 14.8 

-6 04 38 

0.732 

16.73 

16.16 




15.24 

14.72 

14.06 

-26.38 

1321+28 

13 21 15.9 

28 47 19 

0.549 

16.88 

16.15 

16.69 

16.20 


15.39 

14.81 

14.09 

-25.57 

1323+29 

13 23 20.9 

29 10 07 

0.966 

16.91 

16.17 

16.83 

16.70 





-26.84 

1326-05 

13 29 28.6 

-5 31 36 

0.580 

15.47 

15.02 

15.59 



14.51 

13.91 

13.10 

-27.10 

1351+01 

13 51 28.4 

1 03 39 

1.086 

17.68 

16.94 

17.00 

16.50 





-26.31 

1355+02 

13 58 24.0 

2 13 44 

0.955 

16.61 

16.00 







-27.12 

1404+09 

14 04 11.1 

9 37 38 

0.434 

17.22 


17.22 

16.40 





-24.72 

1406-03 

14 06 10.8 

-3 19 13 

0.860 

17.01 

16.21 


16.70 


15.90 

15.63 

14.72 

-26.48 

1415+00 

14 15 49.7 

0 53 56 

1.042 

17.46 

16.91 


16.80 





-26.45 

1415-00 

14 15 28.8 

-0 26 34 

1.148 

17.28 

16.61 


15.90 


15.73 

15.33 

14.61 

-26.84 

1424-00 

14 26 58.5 

-0 20 57 

0.632 

16.52 

15.99 







-26.25 

1528+28 

15 28 40.6 

28 25 30 

0.450 

16.44 

15.79 

16.39 

16.60 


15.42 

14.67 

14.11 

-25.58 

1630+37 

16 30 18.7 

37 19 03 

0.960 

17.00 



17.00 





-26.74 

1633+39 

16 33 02.1 

39 24 28 

1.023 

16.00 

15.50 


15.40 





-27.88 

2241-24 

22 44 40.4 

-24 03 02 

1.958 

16.95 

16.06 

16.83 

16.65 

16.31 

15.95 

15.45 

14.94 

-28.41 

2242-47 

22 45 20.4 

-47 04 50 

0.201 

15.28 

14.22 

14.95 

14.64 

14.13 




-25.07 

2313-30 

23 16 37.9 

-30 29 30 

0.257 

15.62 

15.02 

15.15 

14.97 

14.89 




-25.25 

2321-29 

23 24 26.1 

-28 54 59 

0.279 

16.44 

15.84 

16.03 

15.80 

15.35 

14.75 

13.93 

12.90 

-24.61 









Table 1—Continued 


Name 

R.A.(J2000) 

dec(J2000) 

z 

ms 

mu 

mv 

mR 

mi 

mj 

mn 

niK 

Mb 

2337-33 

23 39 42.7 

-33 31 21 

0.501 

16.91 

16.09 

17.02 

16.66 

16.62 

16.14 

15.42 

14.37 

-25.36 

2352-34 

23 55 25.5 

-33 57 57 

0.702 

16.92 

15.49 

16.45 

16.47 

16.28 




-26.14 

2357-35 

0 00 01.1 

-35 03 35 

0.508 

16.79 

15.48 

16.76 

16.74 

16.51 

16.22 

15.33 

14.51 

-25.51 


'^log W Hz-1 1.4GHz 



Table 2. ISOPHOT Observations LOG 


Source 

Detector 

A 

(nm) 

Int.time 

(s) 

Aperture 

(") 

Chopper 

RASTER 

0034-70 

PI 

7.5 

180 

18 

T 



PI 

11.5 

180 

18 

T 



P2 

25.0 

180 

23 

T 



ClOO 

60.0 

230 

43.5 

T 



ClOO 

100.0 

230 

43.5 

T 


0049-29 

ClOO 

60.0 

720 

43.5 


3x3 


C200 

160.0 

960 

89.4 


2x4 

0051-39 

PI 

11.5 

620 

18 

R 



P3 

60.0 

600 

52 

R 



C200 

160.0 

600 

89.4 

R 


0059-30 

PI 

11.5 

900 

18 

R 



P3 

60.0 

900 

52 

R 



C200 

160.0 

600 

89.4 

R 


0105-26 

PI 

7.5 

230 

18 

T 



PI 

11.5 

230 

18 

T 



P2 

25.0 

230 

23 

T 



ClOO 

60.0 

230 

43.5 

T 



ClOO 

100.0 

230 

43.5 

T 


0120-28 

PI 

11.5 

900 

18 

R 



P3 

60.0 

900 

52 

R 



C200 

160.0 

400 

89.4 

R 


0129-40 

ClOO 

60.0 

720 

43.5 


3x3 


C200 

160.0 

900 

89.4 


2x4 

0144-39 

ClOO 

60.0 

720 

43.5 


3x3 


C200 

160.0 

920 

89.4 


2x4 

0203-08 

PI 

7.5 

330 

18 

T 



PI 

11.5 

330 

18 

T 



P2 

25.0 

330 

23 

T 



ClOO 

60.0 

230 

43.5 

T 



ClOO 

100.0 

230 

43.5 

T 
















Table 2—Continued 


Source 

Detector 

A 

(/im) 

Int.time 

(s) 

Aperture 

(") 

Chopper 

1245-03 

ClOO 

60.0 

520 

43.5 



C200 

160.0 

960 

89.4 


1249-02 

P3 

60.0 

480 

52 

R 


C200 

160.0 

600 

89.4 

R 

1249-06 

P3 

60.0 

330 

52 

R 


C200 

160.0 

360 

89.4 

R 

1252+02 

ClOO 

60.0 

680 

43.5 



C200 

160.0 

860 

89.4 


1321-05 

P3 

60.0 

480 

52 

R 


C200 

160.0 

360 

89.4 

R 

1321+28 

PI 

7.5 

200 

18 

T 


PI 

11.5 

200 

18 

T 


P2 

25.0 

200 

23 

T 


ClOO 

60.0 

230 

43.5 

T 


ClOO 

100.0 

230 

43.5 

T 

1323+29 

PI 

7.5 

230 

18 

T 


PI 

11.5 

230 

18 

T 


P2 

25.0 

230 

23 

T 


ClOO 

60.0 

230 

43.5 

T 


ClOO 

100.0 

230 

43.5 

T 

1326-05 

P3 

60.0 

270 

52 

R 


C200 

160.0 

150 

89.4 

R 

1351+01 

PI 

7.5 

450 

18 

T 


PI 

11.5 

450 

18 

T 


P2 

25.0 

450 

23 

T 


ClOO 

60.0 

230 

43.5 

T 


ClOO 

100.0 

230 

43.5 

T 

1355+02 

P3 

60.0 

330 

52 

R 


C200 

160.0 

350 

89.4 

R 

1404+09 

PI 

7.5 

180 

18 

T 


RASTER 


3x3 

2x4 


3x3 

2x4 










Table 2—Continued 


Source 

Detector 

A 

Int.time 

Aperture 

Chopper RASTER 



(/im) 

(s) 

(") 



PI 

11.5 

180 

18 

T 


P2 

25.0 

180 

23 

T 


ClOO 

60.0 

230 

43.5 

T 


ClOO 

100.0 

230 

43.5 

T 

1406-03 

PI 

7.5 

450 

18 

T 


PI 

11.5 

450 

18 

T 


P2 

25.0 

450 

23 

T 


ClOO 

60.0 

230 

43.5 

T 


ClOO 

100.0 

230 

43.5 

T 

1415+00 

PI 

7.5 

300 

18 

T 


PI 

11.5 

300 

18 

T 


P2 

25.0 

300 

23 

T 


ClOO 

60.0 

230 

43.5 

T 


ClOO 

100.0 

230 

43.5 

T 

1415-00 

PI 

7.5 

300 

18 

T 


PI 

11.5 

300 

18 

T 


P2 

25.0 

300 

23 

T 


ClOO 

60.0 

230 

43.5 

T 


ClOO 

100.0 

230 

43.5 

T 

1424-00 

P3 

60.0 

330 

52 

R 


C200 

160.0 

350 

89.4 

R 

1528+28 

PI 

7.5 

330 

18 

T 


PI 

11.5 

330 

18 

T 


P2 

25.0 

330 

23 

T 


ClOO 

60.0 

230 

43.5 

T 


ClOO 

100.0 

230 

43.5 

T 

1630+37 

PI 

7.5 

200 

18 

T 


PI 

11.5 

200 

18 

T 


P2 

25.0 

200 

23 

T 


ClOO 

60.0 

230 

43.5 

T 









Table 2—Continued 


Source 

Detector 

A 

{fj,m) 

Int.time 

(s) 

Aperture 

(") 

Chopper 

RASTER 


ClOO 

100.0 

230 

43.5 

T 


1633+39 

PI 

7.5 

180 

18 

T 



PI 

11.5 

180 

18 

T 



P2 

25.0 

180 

23 

T 



ClOO 

60.0 

250 

43.5 

T 



ClOO 

100.0 

250 

43.5 

T 


2241-24 

PI 

11.5 

900 

18 

R 



P3 

60.0 

900 

52 

R 



C200 

160.0 

600 

89.4 

R 


2242-47 

ClOO 

60.0 

720 

43.5 


3x3 


C200 

160.0 

600 

89.4 


2x4 

2313-30 

ClOO 

60.0 

900 

43.5 


3x3 


C200 

160.0 

720 

89.4 


2x4 


PI 

11.5 

330 

18 

R 



P3 

60.0 

330 

52 

R 



C200 

160.0 

330 

89.4 

R 


2321-29 

PI 

11.5 

900 

18 

R 



P3 

60.0 

900 

52 

R 



C200 

160.0 

600 

89.4 

R 


2337-33 

PI 

11.5 

900 

18 

R 



P3 

60.0 

900 

52 

R 



C200 

160.0 

600 

89.4 

R 


2352-34 

PI 

11.5 

900 

18 

R 



P3 

60.0 

900 

52 

R 



C200 

160.0 

600 

89.4 

R 


2357-35 

PI 

11.5 

900 

18 

R 



P3 

60.0 

900 

52 

R 



C200 

160.0 

600 

89.4 

R 













Table 3. far-IR Photometry of Quasars" 


Name 

F 7fim 

F 12^m 

F25/im 

F 60/^m 

F lOO/im 

F 160/im 

F 1 RAS 12 

F 1 RAS 25 

F 1 RAS 60 

Firasioo 

0034-70 

<60 

<60 

<90 

<160 

<130 


<60 

<50 

<80 

<210 

0049-29 




<105 


<66 

<60 

<75 

<120 

<300 

0051-39 


<90 


36+18 


438+84 

<120 

<120 

<150 

<450 

0059-30 


<30 


<41 


225+74 

<90 

<120 

<105 

<360 

0105-26 

<100 

<180 

<210 

<100 

<140 


<210 

<120 

170+40 

<300 

0120-28 


<30 


104+60 


<300 

<90 

<90 

180+40 

240+60 

0129-40 




100+30 


<90 

<90 

<120 

<135 

<300 

0144-39 




74+20 


76+26 

<75 

<120 

110+30 

270+70 

0203-08 

<30 

<60 

<100 

<110 

<100 


<60 

<150 

<150 

<400 

1245-03 




87+22 


50+35 

<90 

<240 

<135 

<240 

1249-02 




75+68 


<90 

<120 

<200 

<120 

<300 

1249-06 




217+72 


<150 

<120 

<150 

<100 

<300 

1252+02 




62+25 


76+20 

110+30 

150+50 

<90 

<400 

1321-05 




160+50 


<150 

<100 

<220 

<150 

<220 

1321+28 

<300 

<300 

60+25 

280+130 

290+95 


<90 

120+30 

<110 

<420 

1323+29 

<300 

<90 

<300 

<110 

<150 


<120 

<90 

<150 

<300 

1326-05 




160+80 


<260 

<90 

<210 

100+30 

260+70 

1351+01 

<45 

<90 

60+20 

160+60 

370+130 


<90 

120+30 

<90 

<360 

1355+02 




264+78 


<138 

130+40 

<105 

<80 

200+70 

1404+09 

<90 

<270 

70+25 

380+140 

566+160 


<120 

<105 

<120 

<360 

1406-03 

<45 

<40 

<180 

310+130 

546+110 


<90 

<150 

<160 

<600 

1415+00 

<45 

<60 

54+16 

790+128 

280+80 


<135 

180+40 

<120 

390+100 

1415-00 

<130 

<60 

56+20 

475+120 

190+120 


<99 

<150 

140+40 

280+100 

1424-00 




69+35 


<120 

<120 

<120 

<135 

<210 

1528+28 

<60 

<30 

50+40 

180+140 

<120 


<80 

<60 

<120 

<220 

1630+37 

<45 

<54 

63+20 

<180 

150+120 


<60 

< 70 

< 90 

<270 

1633+39 

<200 

<240 

<120 

<100 

<100 


<60 

<75 

< 75 

<120 

2241-24 


<24 


45+23 


141+83 

< 90 

<150 

<90 

<240 

2242-47 




<62 


<40 

<90 

<60 

<120 

<360 

2313-30 


<30 


77+20 


<50 

90+20 

<150 

130+30 

<600 


98±84 <100 



Table 3—Continued 


Name 

F 12 /im F 25/im 

F 60/^m F lOO/^m 

F 160/im 

F 1 RAS 12 

F 1 RAS 25 

F 1 RAS 60 

Firasioo 

2321-29 

<20 

110±50 

280±80 





2337-33 

<27 

<36 

<80 





2352-34 

<18 

<33 

<76 

<90 

<90 

<100 

<300 

2357-35 

<30 

<30 

290±72 

140T30 

<120 

140T30 

230T70 


^Fluxes are in mJy 







